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Introductory  rgnarkg 

The  contract  was  awarded  in  1964  and  provided  financial  support  for 
three  years.  Vtark  began  in  May  1964  and,  in  the  subsequent  X  months,  two 
rigs  have  been  constructed  and  used  for  the  experimental  inveetigations 
described  in  references  1  to  5  and  a  third  rig  constructed  in  preparation 
for  future  Mark.  The  three  rigs  simulate  aspects  of  the  Haws  in  gin 
barrels  and  represent  a  progression  towards  gins  so  that,  for  example,  the 
second  rig  removed  problems  associated  with  the  measurement  of  projectile 
velocity  and  the  third  rig  will  increase  the  maximait  projectile  velocity 
to  500  m/s.  The  research  is  part  of  a  larger  program  involving  staff  of 
Imperial  College,  the  Fraunhofer- Institute  fuer  Kurzeitdynamik  (EMI-AFB) 
and  the  Ballistic  Research  laboratory  (BHL)  at  Aberdeen  Proving  Ground. 
Close  contact  with  BU-AFB  and  HRL  has  been  maintained  throughout  the 
contract. 

The  measurements  reported  and  referenced  were  obtained  to  improve 
understanding  of  the  physical  phenomena  associated  with  gun  barrel 
configurations  but  without  the  complications  of  combustion  and  heat 
transfer.  Thus  emphasis  was  placed  cn  fluid  dynamic  properties  including 
pressure  and  velocity  characteristics  in  single-  and  two-phase  unsteady 
flows.  The  results  are  presented  in  sufficient  detail  to  aid  the 
development  and  appraisal  of  calculation  methods  which  involve  the 
numerical  solution  of  conservation  equations  in  finite-difference  form. 
Since  the  calculation  of  flows  involving  two  phases,  requires  that 
corresponding  models  be  devised,  experiments  were  carried  out  in  simpler 
but  relevant  flow  configurations  and  are  described  in  reference  6. 
Further  experiments  were  performed,  partly  supported  from  the  sources,  to 
inprove  techniques  for  the  maasuranent  of  two- phase  flows  (7  to  10). 

A  considerable  part  of  the  overall  effort  was  devoted  to  the  design 
and  testing  of  the  three  rigs  and  of  relatsd  instrumentation .  The 
progression  towards  the  third  rig  involved  the  development  of  mechanists 
and  procedures  for  securing,  releasing,  catching  and  relocating  the 
projectile,  methods  for  measuring  the  projectile  velocity  ever  a  wide 
range  of  values  and  a  configuration  which  would  allow  the  introduction  of 
a  solid  phase  in  particle  form.  In  addition,  laser  velocimetry  was 
developed  to  measure  the  velocity  characteristics  of  particles  ever  an 
extansive  range  of  values  and  with  accuracy  appropriate  to  the  dual 
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purposes  of  inproved  understanding  of  the  flow  and  of  calculation  methods. 
These  aspects  of  the  work  are  not  considered  in  this  report  and  it  should 
be  emphasised  that  they  were  essential  to  the  determination  of  the  results 
and  provide  the  basis  for  the  intended  future  work  with  the  high-speed 
rig. 


The  results  obtained  in  the  three  gin-simulation  rigs  are  discussed 
briefly  in  the  following  section  and  further  details  are  available  in  the 
references  cited.  The  report  ends  with  a  sarmary  of  the  more  important 
conclusions. 

Results  obtained  in  gun- simulation  experiments 

Preliminary  measurements  in  the  first  rig,  see  figure  1,  indicated  the 
main  features  of  in-cylinder  flows  for  initial  pressures  up  to  9  bars. 
The  near-wall  measuronents  revealed  thin  boundary  layers  and  the  maximan 
velocity  fluctuations  in  the  bulk  flow  grew  frcm  zero  to  around  6%  as 
shown  in  figure  2.  They  were  the  first  available  measurements  for 
in-cylinder  flaws  caused  by  a  moving  projectile  and  were  used  at  BKL  to 
support  the  development  of  the  t  code.  They  also  showed  that  improvements 
could  be  made  to  the  rig  and  the  instrumentation  so  as  to  permit  more 
accurate  results  at  higher  pressures.  The  instrumentation  used  with  the 
first  rig  included  a  pressure  transducer,  laser  velocimeter  and  optical 
encoder,  but  the  transducer  signal  was  not  digitised  or  multiplexed,  the 
velocimeter  was  not  arranged  specifically  for  the  gun  simulator  and  the 
encoder  operated  frcm  a  mechanical  rack  and  pinion. 

The  second  rig  was  designed  and  constructed  to  allow  more  detailed 
measurements  in  the  near-wall  region  and  operation  without  the  drag  force 
associated  with  the  rack-and-gear  mechanism.  Improvements  were  also  made 
to  involve  higher  pressures,  to  allow  the  introduction  of  particles  in  a 
quantitatively  known  way  and  to  allow  the  release  of  the  projectile  in  a 
manner  which  could  also  be  used  in  the  third  rig,  proposed  for  use  with 
much  higher  pressures.  The  method  of  measuring  the  projectile  velocity 
was  changed  to  make  use  of  an  optical  technique  and  the  laser-velocimeter 
was  rearranged  with  higher  power  laser  and  faster  signal  processing;  these 
two  electro-optic  devices  wmre  designed  so  that  they  could  also  be  used 
with  the  third  rig.  The  second  rig  is  shown  in  figure  3. 


3 


The  measurements  of  reference  3  and  4  were  obtained  in  the  second  rig 
and,  as  can  be  seen  from  figure  4  for  example,  the  near-wall  velocity 
traces  exhibit  increasing  velocity  fluctuations  as  the  wall  is  approached . 
The  ensemble-averaged  boundary- layer  profiles  follow  clearly  a  power  law 
with  index  varying  from  6.6  to  9.9  as  the  projectile  progresses  dowi  the 
barrel .  These  results,  and  the  calculated  shape  factors  tAuch  varied  frcvn 
1.2  to  1.3,  confirm  that  the  boundary- layer  is  turbulent  with 
characteristics  similar  to  those  of  a  steady  turbulent  boundary  layer. 
The  results  of  reference  3  provide  a  well-defined  picture  of  the  flow  in 
the  simulated  gun  with  an  initial  condition  corresponding  to  quiescent 
flow.  They  were  extended,  in  reference  4,  to  involve  single-  and 
two-phase  flews,  which  involve  flow  into  the  breech  through  an  array  of 
jets,  and  confirm  that  the  geometry  of  the  initial  chamber  does  not  affect 
the  downstream  flow  properties. 

The  two-phase  flow  experiments  involved  500  ym  diameter  particles  with 
initial  volume  concentrations  up  to  1.3%.  A  sample  of  the  results  is 
shown  on  figure  5  and  allows  comparison  with  those  obtained  without  the 
particles.  The  results  were  obtained  at  the  first,  second  and  third 
windows  corresponding  to  axial  locations  of  4.%,  7.90  and  10.84d.  It  can 
be  seen  that  the  solid  particles  lag  behind  the  fluid  by  times  which 
increase  with  axial  distance,  i.e.,  as  the  projectile  is  approached.  The 
difference  between  fluid  and  particle  velocity  increases  in  the  bulk  flew 
from  about  1  m/s  and  x/d  =  4.96  and  t  =  20.5  ms  to  8  m/s  (21%)  at 
x/d  =  10.84  at  t  =  27.5  ms.  The  lag  is  due  to  inertial  effects  associated 
with  solid  particles  in  a  strongly  accelerated  flow  field  and  can  be 
expected  to  be  much  greater  at  the  accelerations  of  the  third  rig  and  of 
guns. 


The  work  of  reference  5  was  carried  out  to  determine  the  effect  of  a 
three- fold  increase  in  particle  concentration,  to  re-examine  the  near-wall 
region  and  to  test  procedures  vAiich  are  to  be  used  in  the  third  rig.  The 
measurements  in  the  near-wall  region  have  confirmed  the  turbulent  nature 
of  the  wall  boundary  layer  and  the  increase  in  particle  concentration  has 
shown  a  tendency  for  the  particles  cn  the  centre  line  to  decrease  in 
velocity  less  rapidly  than  with  the  lower  concentrations  examined 
previously.  The  instrunentation  and  experimental  techniques  proved 
satisfactory  and  led  to  the  arrangonent  of  the  third  rig,  as  shown  on 
figure  6.  Preliminary  tests  with  this  rig  have  confirmed  that  the 
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projectile  velocity  corresponds  to  values  determined  with  isentropic-flow 
assumptions ,  as  shown  on  figure  7,  and  that  values  of  over  400  m/s  can 
readily  be  attained. 

The  preparations  fear  future  work  include  the  provision  of  a  high-speed 
counter,  an  argon-ion  laser  and  a  Bragg  cell  optical  unit.  The  software 
has  been  prepared  and  provision  made  for  the  necessary  interface  between 
the  centre  and  IBM  AT  computer.  In  addition,  the  related  efforts  of 
re''*-.-'  _e  6  has  quantified  the  influence  of  particles  in  liquid  with 
volume  concentrations  up  to  around  0.15  and  references  7  to  10  provided 
instrumentation  for  the  simultaneous  maasuronent  of  the  discrete  and 
carrier  phase  with  known  accuracy.  Work  is  in  progress  to  determine  the 
accuracy  with  vdiich  the  barrel  temperature  can  be  measured  as  the 
projectile  passes  downstream. 

Conclusions 


The  following  paragraphs  provide  a  summary  of  the  more  important 
conclusions: 

1.  The  first  and  second  gun  simulators,  and  related  instrumentation,  have 
provided  information  of  the  velocity  characteristics  of  projectile  and 
gas  as  a  function  of  projectile  position  and  initial  pressure .  The 
experiments  have  provided  guidance  for  the  design  arid  construction  of 
a  third  rig,  and  related  instrumentation,  which  will  allow  projectile 
velocities  up  to  500  m/s. 

2.  The  relationship  between  breech  pressure  and  projectile  velocity  is 
shown  to  be  described  by  one-dimensional,  unsteady  forms  of  the 
continuity  and  momentum  equations .  The  measured  record  of  braech 
pressure  leads  to  calculated  values  of  projectile  velocity  within  0.5% 
of  the  measured  values  and  confirms  that  friction  between  projectile 
and  tube  well  is  negligible. 

The  velocity  results,  for  single-  and  two-phase  flow®,  are  provided  in 
sufficient  detail  to  aid  the  development  and  appraisal  of  calculation 
methods,  based  on  the  numerical  solution  of  conservation  equations  in 
differential  form.  They  quantify  the  variation  of  gas  and  particle 
velocity  as  a  function  of  position  and  time  and  quantify  its  decay 
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after  the  passage  of  the  projectile.  The  well  boundary  layer  is  shown 
to  be  turbulent  and  the  lag  between  500  pm  solid  particles  and  the  gas 
phase  quantified  for  two  values  of  concentration . 

4.  The  effect  of  initial  charter  geometry  is  shown  to  be  snail  at  the 
aamapratively  low  velocities  of  the  two  simulators. 

5.  The  experimental  and  measurement  techniques  have  been  shown  to  be 
appropriate  and  will  be  used  in  a  subsequent  investigation  in  a  third 
rig  with  much  higher  pressures  and  velocities.  Ibis  rig  has  been 
constructed  and  tested.  Preliminary  measurements  confirm  conclusion  2 
and  show  that  it  will  be  possible  to  extend  the  research  to  projectile 
velocities  in  excess  of  500  m/s,  which  overlaps  with  the  lower  limit 
of  the  related  experimental  program  at  the  Ernst  Mach  Institute.  This 
third  rig  has  been  designed  to  permit  measurements  with  pellets  to 
solid  propellant  to  assist  the  acceleration  of  the  projectile  and 
related  two-phase  flow  instrumentation  has  been  developed  and  tested . 
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Figure  1 


Arrangement  of  first  projectile  rig 
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Figure  2  Centre-line  velocity  traces 
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Figure  3 


Arrangement  of  second  projectile  rig 
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Figure  4 


Velocity  variation  with  time  and  distance  from  the  wall 
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Appendix 

The  imegral  forms  or  the  momentum  and  energy  solutions 
arc  obtained  below  for  the  case  where  freestream  vorticity 
and  total  enthalpy  gradients  are  included  in  the  boundary- 
layer  solution.  From  Eqs.  (I)  and  (2).  and  (8). 

J  |W7(  ifiU)m-pu)  4 

*i Cl(pv)m,-0v}  4  | dr  (Al) 

Equation  (Al)  together  with  Eqs.  (I>,  00),  and  Ha)  then 
gives 

r„«  (u4 f  [  (pp}"-pa]di 

dr*  dr* 

-  — -u),',  [(pu)m -pu]yd\  +  — -l  pu(u"‘  -u)dy 

dv  Jm  dr  -’«• 

(A2) 

Similarly,  from  Eqs  (I).  (4).  (7b),  and  (10), 
q»-  [(  [H"'{  {pu)n,-pu}  [pu(H"  - H) J, 

4  H:\ (pc)"'  -pel  4  \pv(H"  -//)),  |d.» 

f  *  d 

*  ( ( pu ) —  pu  J d„v  4  —  ( H  );. 

.««  dr 

f A  d  f  * 

x  |(p</)"'-pu).vd»  4  — pu{H  -H) d.t  (A3) 
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Introduction 


VELOCITY  and  pressure  have  been  measured  in  the  flow 
behind  a  projectile  traveling  in  a  tube  and  propelled  by 
compressed  gas  The  arrangement  is  intended  to  simulate  that 
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of  interior  ballistics,  with  non  reacting  flow  and  projectile 
velocities  up  to  21  m/s.  Although  the  real  processes  are  much 
more  complex  and  involve  a  continuous  interaction  of  two 
phases  (solid  gas  or  liquid  gas)  through  a  combustion  process, 
the  present  simplified  experiments  serve  to  provide  fundamen¬ 
tal  understanding  of  the  flow  behind  in-bore  projectiles  and 
thereby  support  the  development  of  phenomenological 
models,  such  as  those  of  Refs.  I  and  2,  and  multidimensional 
solution  methods  such  as  that  of  Ref.  3. 

A  projectile  was  secured  in  position  inside  a  tube  pressur¬ 
ized  by  compressed  gas  to  a  required  pressure.  After  the  gas  had 
become  quiescent,  the  projectile  was  released.  The  pressure  in 
the  lube  and  the  projectile  velocity  were  measured,  together 
with  the  local  velocity  of  the  fluid  obtained  with  a  laser 
velocimeter,  as  a  function  of  lime.  The  experiment  was  re¬ 
pealed  at  different  locations  until  a  satisfactory  picture  of  the 
flow  had  been  assembled.  The  initial  pressure  was  regarded 
as  a  variable  so  as  so  determine  the  extent  to  which  the  flow 
properties  varied  with  projectile  velocity,  and  the  experiment 
was  conducted  with  (wo  initial  chambers  of  different  length 

Experimental  System 

The  flow  arrangement  and  related  instrumentation  are 
shown  in  Fig.  I .  The  measurements  were  obtained  tn  the  cavity 
formed  by  a  76.7-mm-diam  tube,  a  blanked  end.  and  a 
76  4-mm-diam  projectile  (hat  had  a  flat  end  The  effect  of  ini¬ 
tial  chamber  length  was  investigated  by  using  two  different 
lengths  of  177.3  and  31 1  mm,  which  were  achieved  by  altering 
the  projectile  overall  length.  The  upper  300  mm  of  the  tube 
was  made  of  plexiglass  and  the  remaining  pan  of  mild  steel 
with  the  gap  between  the  projectile  and  tube  wall  sealed  by  two 
silicon -rubber  rings.  The  projectile  was  allowed  to  travel  for  a 
distance  of  approximately  400  mm  before  being  retarded  by  a 
series  of  compression  springs. 

The  initial  gas  pressure  P,  ranged  from  2.91  to  9.1  bars  and 
was  achieved  with  nitrogen  taken  from  a  pressurized  cylinder 


Fig.  I  Schematic  diagram  of  experimental  omem. 
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Fig.  2  CtMffliM  fflocH)  tram  at  z-  153.4  nun  for  two  fait t* I 
cbaraber  leogthi  awl  wriwt  Mtial  pttium 


t  ( a»(  I 

Fig.  3  CMrHItr  velocity  tram  at  various  axial  locatfoat  for  ikon 
kohial  chamber  a  ad  P,  •  4.2  ban. 


The  initial  pressure  and  the  variation  of  pressure  with  time 
were  measured,  with  a  precision  of  ±.002  bar.  using  a 
piezoelectric  transducer  and  a  charge  amplifier  whose  output 
was  digitized  and  interfaced  to  a  microprocessor 

Laser  Doppler  anemometry  was  used  to  obtain  the  fiu;d 
velocities  with  the  optica)  arrangement  of  Ref  4.  Silicon  oil 
droplets,  served  as  light -scattering  particles,  were  introduced 
in  the  initial  chamber  prior  to  the  charging  of  the  volume  w  ith 
(he  pressurized  gas 

The  frequency  of  the  Doppler  signal  was  measured  by  a 
purpose-built  frequency  counter  also  interfaced  to  the 
microprocessor.  The  time  information  corresponding  to  each 
frequency  count  was  obtained  from  a  Programmable  Time 
Module  initiated  to  time  at  a  rate  of  4P  KHz  after  the  fir* 

I -mm  travel  of  the  projectile 

The  projectile  position  was  determined  by  an  optical  en¬ 
coder  coupled  to  the  rack  and  gear  arrangement  of  Fig.  1 .  and 
each  encoder  pulse  corresponded  to  an  axial  distance  of 
0.2  mm.  The  lime  between  two  consecutive  pulses  was 
measured  by  a  high-frequency  crystal  clock  oscillating  at  a  rate 
of  1.0  MHz,  and  values  were  stored  in  a  buffer  memory. 

The  data  stored  were  subsequently  processed  by  the  micro¬ 
processor  to  determine  and  display  the  variation  of  projectile 
and  local  flow  velocity  and  pressure  as  a  function  of  time.  The 
accuracy  of  fluid  velocity  measurement  is  expected  to  be  of  the 
order  of  1%  except  in  the  near-wall  region  where  velocity - 
gradient  broadening  became  important  and.  with  the  present 
optical  arrangement,  precluded  useful  results  within  around 
0.5  mm  of  (he  wall. 

Rests**  and  Discussion 

The  traces  of  center  hue  velocity  Ut  measured  at  an  axial 
distance  of  153.4  mm  from  the  closed  er  d,  are  shown  in  Fig  2 
for  both  chambers  and  for  a  range  of  initial  pressures.  The 
velocities  with  the  larger  initial  volume  are  lower,  since  the 
•Kasurement  location  is  farther  from  the  projectile.  It  is  also 
evident  that  each  velocity  trace  is  initially  associated  with  a 
monotonic  increase  and  that  fluctuations  gradually  appear 
with  decreasing  acceleration.  The  root-mean-square  value  of 
these  fluctuations  is  around  3%,  with  normalized  peak-to- 
peak  values  of  about  6%. 

Centerline  velocity  traces  are  also  shown  in  Fig.  3  for  the 
short  initial  chamber  and  initial  pressure  of  6.2  bars  The 


t  iols  of  ooodimtMioMl  cvwieriloc  velocity  f  ogoiosl  ooo- 
'Ml  BXiol  dhtance  Z*  for  long  chamber  a*4  P,  •  2.41  bore: 
c/b*.  or  -  It.  *r»)t.  or- M.  bW  » t »  40  mt:  for  long  chamber 
ooi  P,  •  4.2  bon;  or- 5,  ■i-IO.mM  * r -  2t  rat:  for  *on  chamber 
mi  Pt*M  bon;  •  /- J.  or-  It  mi  «/•»■>. 

results  include  some  of  those  of  Fig.  2.  extend  to  larger  values 
of  axial  distance,  and  include  positions  initially  occupied  by 
the  projectile.  The  fluid  velocity  increases  with  axial  distance 
and  in  accord  with  the  accelerating  projectile  At  positions  ini¬ 
tially  occupied  by  the  projectile,  the  velocity  decays  after  the 
projectile  has  passed,  ai  a  rate  that  increases  with  the  axial 
distance.  The  results  confirm  the  tendency  of  the  fluid 
velocities  to  fluctuate  with  an  amplitude  (hat  appears  to  in¬ 
crease  until  a  peak-to-peak  value  of  around  6%  is  achieved 

Figure  4  shows  the  nondimensionat  centerline  velocity 
C P  ( (J,  /(J- 1  plotted  against  the  nondimensiona!  axial  distance 
*•  U%.  I  for  different  initial  pressures,  chamber  lengths,  and 
times  to  the  cycle;  U?  and  i0  are  the  corresponding  projectile 
velocity  and  position,  respectively.  The  figure  clearly 
demonstrates  that  the  fluid  velocity  distribution  between  the 
breech  and  the  projectile  bate  is  linear  and  scales  with  the  pro¬ 
jectile  velocity  at  all  times  in  the  cycle,  for  the  range  of  initial 
pressures  and  two  chamber  lengths  considered  here. 

Radial  profiles  of  velocity  were  obtained  from  velocity 
traces  such  as  those  of  Fig.  5  and  were  found  to  be  nearly 
uniform  apart  from  thin  regions  near  the  wall  boundaries.4 
The  traces  of  Fig.  5  correspond  to  locations  0.5, 1,  and  2  mm 
from  the  wall  with  the  velocity  values  at  rm  ■  0.5  mm  generally 
of  the  order  of  90-94%  of  the  velocities  at  r  „ «  2  mm.  The 
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Fig-  5  Fluid  retoell.t  men  nt  wall  boulder.  «:•  15J.4  mm  lot 
6«n  lolltel  chamber  ood  P,  -  6.2  ban. 
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boundary-layer  thickness  was  found  lo  be  less  than  1  mm. 
which  corresponds  to  about  I of  the  tube  diameter.  The 
numerical  calculations  of  Schmitt  et  el.’  predict  similar 
boundary-layer  thickness.  The  turbulence  levels  are  slightly 
lower  than  those  observed  in  the  core  region. 

Conclusions 

The  main  findings  of  this  investigation  are  as  follows: 

1)  The  fiow  at  any  given  time  in  the  cycle  is  nearly  one- 
dimensional  apart  from  thin  regions  near  the  wall  boundaries. 
Velocity  fluctuations  reach  a  maximum  peak-to-peak  value  of 
around  6 We  in  the  core  region. 

2)  In  all  cases  the  boundary-layer  thickness  is  less  than  I  mm. 
which  corresponds  to  about  1%  of  the  bore.  The  velocity  fluc¬ 
tuations  in  the  boundary  layer  are  generally  lower  than  those  in 
the  core  region. 

3)  The  velocity  increases  linearly  with  axial  location  from 
aero  at  the  breech  to  that  of  the  projectile  at  the  base. 

4)  The  bulk  fluid  velocity  scales  with  that  of  the  projectile,  in¬ 
dependently  cf  'i.itial  chamber  length  and  pressure. 
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abstract 

Ad  experimental  investigation  of  the 
fluid  motion  in  the  wake  of  an  in¬ 
cylinder  projectile  accelerated  from 
rest  by  compressed  gas  is  reported. 
Time-resolved  measurements  of  velocity 
were  obtained  by  laser  Doppler  anemo- 
metry  as  a  function  of  time  and  for  a 
range  of  positions  in  the  flow. 

The  assembly,  wbicb  simulates  a  gun 
barrel,  comprised  a  plexiglass 
cylinder  and  a  flat-based  projectile  of 
76.4  mm  diameter.  The  initial  gas 
pressures  ranged  from  1.7  to  3.6  bars 
and  resulted  in  maximum  projectile 
velocities,  over  a  300  mm  travel,  from 
4.5  ms'1  to  9.2  ms*  .  The  correspond¬ 
ing  Reynolds  number  range,  based  on 
cylinder  diameter  and  atmospheric 
density,  was  2.3  to  4.7  104  . 

The  results  indicate  that  the  flow  is 
near  one-dimensional  with  boundary 
layers  generally  of  the  order  of  17 
of  the  cylinder  diameter  and  maximum 
turbulence  intensities  of  about  47. 
Velocity  traces  obtained  near  tbe  wall 
boundary  do  not  indicate  turbulent 
fluctuations  suggesting  that  the 
boundary  layers  remain  laminar.  The 
results  also  reveal  that  the  spatial 
distribution  of  velocity  between  the 
breech  and  projectile  base  is  linear 
at  any  given  time  in  tbe  cycle. 


NOMENCLATURE 

D  Cylinder  bore 

Pj  Initial  chamber  pressure 

ra  Radial  distance  from  cylinder  wall 

t  Time 

D  Axial  velocity  of  fluid 

Dg  Centreline  velocity 

Up  Projectile  velocity 

UpB4x  Maximum  projectile  velocity 
u  Rats  of  axial  ve'.oelty  fluctuations 


u/U  Turbulence  intensity 

2 

Vj  Initial  chamber  volume  («D  Zj/4) 
Initial  chamber  length 
z,r  Cylindrical  coordinates 
zp  Projectile  position  (Zj  *  zt ) 
zt  Projectile  travel 
i  Boundary  layer  thickness 
£d  Displacement  thickness 

u  Viscosity  of  fluid 
p  Density  of  fluid 


1  INTRODUCTION 


The  processes  which  take  place  in 
interior  ballistics  are  complex  and 
involve  a  continuous  interaction  of  twc 
phases  (solid-gas  or  liquid-gas)  through 
a  combustion  process  which  converts  the 
chemical  energy  into  heat.  A  number  of 
phenomenological  models  have  been 
developed,  see  for  example  reference  1 
and  2,  to  predict  these  processes  and 
subsequent  projectile  motion.  They  have, 
however,  bees  inadequately  tested  and  the 
need  to  solve  problems  such  ss  tube  wear 
requires  a  more  detailed  and  assured 
understanding  of  the  physical  processes 
taking  place,  particularly  close  to  the 
tube  walls.  Tbe  representation  of  these 
processes  requires  a  multi-dimensional 
approach  such  ss  that  of  reference  3 
which,  is  turn,  requires  that  tbe 
physical  assumptions  be  confirmed  and 
quantified  by  experiments. 

The  present  work  is  designed  to  provide 
the  experiments!  information  in  the  form 
of  fluid  velocity  characteristics  for  a 
single-phase  flow  and  is  the  absence  of 
combustion.  Although  tbe  velocities 
involved  were  an  order  of  magnitude 
smaller  than  those  found  la  practice, 
the  work  provides  a  basic  un  lerstasding 
of  the  ls-cyllnder  flow  behind  a 
projectile,  directly  relevant  to  interior 


2 


ballistics. 

In  tbs  following  section  the  test  rig 
Is  described  together  with  the  experi¬ 
mental  techniques  used  for  projectile 
sod  fluid  velocity  measurements  and  a 
brief  discussion  on  the  experimental 
uncertainties.  Results  are  presented  and 
dlacuaaed  in  Section  3, and  Section  4 
summarises  the  more  important  findings. 


2  EXPERIMENTAL  STBTEH 

2.1  Flow  Configuration 

A  achematlc  diagram  of  the  experimental 
system  is  shown  in  Figure  1 .  The 
assembly  was  mounted  vertically  on  a 
rigid  table  and  comprised  a  flat -based 
projectile  of  76.4  mm  diameter,  driven 
inside  a  mild  steel  and  plexiglass 
cylinder  of  76.7  mm  diameter  a6  indicated 
in  Figure  1.  The  main  regions  of  interest 
included  the  initial  chamber  and  part  of 
the  volume  swept  by  the  projectile,  and 
corresponded  to  the  plexiglass  cylinder 
which  provided  the  optical  access  for  the 
laser  Doppler  anemometer.  The  plexi¬ 
glass  cylinder  was  296.3  mm  in  length  and 
had  a  wall  thickness  of  12.3  mm.  The 
clearance  gap  of  0.15  am  between  the 
cylinder  wall  and  projectile  was  sealed 
in  the  metal  part  of  the  cylinder  by 
means  of  two  silicon  rubber  rings  as 
shown  in  Figure  1.  The  friction  between 
the  cylinder  and  projectile  was  kept  to  a 
minimum  by  PTFE  guide  rings  lubricated  by 
graphited  grease.  The  total  travel 
distance  available  before  the  projectile 
comes  to  rest  was  about  400  am  aod  a 
system  of  compression  springs  retarded 
the  projectile  before  it  ban  travelled 
this  distance;  this  was  necessary  in 
order  to  reduce  the  Impact  loading  on  the 
projectile  and  end  plate.  The  travel 
distance  unaffected  by  the  spring 
retardation  was  approximately  300  am.  The 
main  geometric  characteristics  of  the 
assembly  are  stumsarised  in  Table  1. 

TABLE  1  GEOMETRIC  CHARACTERISTICS  OF  THE 
GUN  ASSEMBLY 


cylinder  by  means  of  initially  compressed 
nitrogen  gas  at  ambient  temperature.  The 
bottled  nitrogen  was  chosen  as  the  work¬ 
ing  fluid  as  it  was  found  to  be  dryer 
than  other  coaaerclally  available 
bottled  gases  and  as  humidity,  at  temper¬ 
atures  below  -B0°C,  causes  condensation 
with  the  resulting  water  droplets 
contributing  as  noise-induced  counts  in 
the  laser  Doppler  anemometry  measurements 
Initial  gas  pressures  P,  were  in  the 
range  of  1.7  to  3.B  bars  and  resulted, 
over  the  300  mn  travel ,  maximum 
project llejvelocities  fD  ,  1  of  4.5  ms"3 
to  9.2  ms'3.  The  corresponding  Reynolds 
number  range,  based  on  the  cylinder 
diameter  and  atmospheric  density 
(D1 _ pP/u),  was  2.3  to  4.7  10*. 

2.2  Data  Acquisition  and  Processing 

An  optical  rotary  encoder  coupled  to  a 
gear  was  used  is  conjunction  with  a 
rack  mechanism  driven  linearly  by  the 
projectile  to  determine  its  dynamic 
characteristics.  The  encoder  provided 
1000  pulses  per  revolution  and  with  the 
particular  gear  diameter  used  (63  mm) 
resulted  in  a  linear  displacement  of 
0.2  nan  per  pulse.  The  time  between  two 
consecutive  pulses  were  determined  by 
a  high  frequency  crystal  clock  oscillat¬ 
ing  at  a  rate  of  1.0  MHz  and  values  were 
stored  is  a  buffer  memory  for  subsequent 
processing  in  a  microprocessor  to  yield 
the  velocity-travel  characteristics  of 
the  projectile. 

Laser  Doppler  anemometry  was  used  to 
obtain  the  fluid  velocities  with  an 
optical  arrangement  similar  to  that  used 
in  reference  4.  The  anemometer  made  use 
of  forward-scattered  light  and  comprised 
s  low  power  Be-Ne  laser,  a  focussing 
lens  L.,  a  rotating  diffraction  grating 
which  Splits  the  beam  and  provides  the 
frequency  shift,  a  collimating  lena  Lj, 
another  focussing  less  L* ,  collecting 
less  L-  aod  a  photomultiplier  which 
converts  the  optical  signal  into 
electrical  signal.  The  principal 
characteristics  of  the  optical  system 
used  are  gives  is  Table  2. 


Bore  diameter,  D 
Initial  chamber  length,  Z, 
Initial  chamber  volume,  V£ 
Maximum  travel  of  projectile 
Unaffected  travel  of 
projectile 


76.7  so 
177.3  mm 
0.82  ltr 
400  me 

300  bd 


The  projectile  was  impulsively  started 
from  rest  using  a  nlnlmum-force-bearlqg 
release  mechanism  as  indicated  is 
Figure  1,  and  accelerated  inside  the 


TABLE  2  PRINCIPAL  CHARACTERISTICS  OF 
THE  OPTICAL  SYSTEM 


Half  angle  of  beam 
intersect  ion 

Fringe  spacing 

Length  of  probe  volume 
at  e"*  intensity  level 

Diameter  of  probe  volume 
at  e*z  lntensl.y  level 


6.1  degrees 
2.98  im> 


0.58  mm 


82 


3 


Number  of  fringes  within 
t*  intensity  level  21 

Frequency  shift  1.0  KHz 


Silicon  oil  particles  were  employed  ss 
light  scattering  particles  and  were 
obtained  by  using  low  pressure  nitrogen 
is  a  blast  type  atomiser  as  described  in 
reference  4.  Particles  were  introduced 
in  the  initial  chamber  prior  to  the 
charging  of  the  volume  with  the  pressur¬ 
ised  gas.  The  initial  particle 
concentration  obtained  by  this  method  was 
sufficient  to  obtain  the  measurements  of 
time-resolved  velocities  at  almost  every 
part  of  the  cylinder. 

The  Doppler  signals  from  the  photo¬ 
multiplier  output  were  band-pass 
filtered  (0.2  -  5.5  KHz  at  -3  dB)  to 
remove  the  pedestal  and  the  high 
frequency  noise  prior  to  their  input  to 
a  purpose-built  frequency  counter  with  a 
maximum  sampling  rate  of  about  22  kHz. 
although  there  was  no  direction 
ambiguity  in  the  fluid  velocities 
encountered,  a  minimum  frequency  shift 
of  1.0  MHz  was  used  to  ensure  that  more 
than  32  fringes  were  present  in  each 
Doppler  burst  as  required  for  efficient 
frequency  counting,  see  reference  5.  The 
input  and  output  signals  of  the  counter 
were  input  to  an  oscilloscope,  as  shown 
in  Figure  1,  to  monitor  the  signal 
quality  and  to  ensure  its  satiafactcr; 
operation  during  each  shot .  Th« 
frequency  counter  was  interfaced  to  tbt 
microprocessor  which  sampled  the  digital 
output  of  the  counter  and  stored  it  in  the 
memory  via  a  6502  machine  code  at  a 
maximum  rate  of  17  kHz.  The  time 
information  corresponding  to  each 
frequency  count  was  obtained  from  a 
Programmable  Timer  Module  (PTHJ  which  was 
initiated  to  time  at  a  rate  of  417  kHz 
with  the  arrival  of  the  5th  encoder 
pulse,  corresponding  to  the  first  1  mm 
displacement  of  the  projectile. 

The  data  stored  lo  the  microprocessor  and 
buffer  memories  were  then  processed  by  a 
software  written  in  "BASIC"  to  determine 
both  the  projectile  and  fluid  dynamics  in 
aacb  shot,  and  subsequently  to  plot  and 
display  them  on  a  VDU. 


2.3  Precision 

For  the  measurements  of  projectile 
velocity  a  clock  frequency  of  1.0  MHz  was 
used ,  as  mentioned  above,  to  time  the 
projectile  motion  during  each  encoder 
pulse  corresponding  to  an  axial  distance 
of  0.2  an.  Due  to  the  count  ambiguity 
associated  with  this  clock  frequency, 


there  was  an  Inherent  error  in  the 
meaauremed  projectile  velocities,  which 
increased  with  increasing  velocity.  The 
maximum  error  In  the  velocity  range 
encountered  was  about  3%,  but  the  overall 
error  is  expected  to  be  much  less  in 
most  parts  of  the  cycle. 

The  resolution  of  the  frequency  counter 
is  frequency  dependent  and  for  the 
frequency  range  encountered  (1-3  MHz) 
corresponds,  on  average,  to  4.7  kHz  , 
which  Implies  a  resolution  of  0.014  ms' 
in  terms  of  fluid  velocity.  The 
associated  error  for  as  average  fluid 
velocity  of  3  ms”1  was,  therefore,  less 
than  0.51.  The  rms  of  the  short  term 
instability  of  the  frequency  shift  was 
about  0.2 T  which  implies,  together  with 
the  resolution  of  the  frequency  counter, 
an  overall  error  of  no  more  than  IT  in 
the  measured  velocity  values. 

The  initial  chamber  pressure  was  read  os 
a  pressure  gauge.  Preliminary  checks 
made  by  mercury  manometers  up  to  1.5  bars 
indicated  maximum  deviations  of  less 
than  f  21  from  the  readings  of  the 
pressure  gauge. 

These  magnitudes  of  the  uncertainties 
considered  above  are  not  expected  to 
affect  the  conclusions  reached  on  the 
basis  of  the  results  which  are  presented 
and  discussed  in  the  following  section. 


3  RESULTS  AMD  DISCDSSION 

The  results  presented  and  discussed  in 
this  section  are  representative  samples 
chosen  to  describe  the  flow  process. 
Although  some  results  are  presented  for 
maximum.projectlle  velocities  from 

4.5  ms”1  to  B.2  ms”1  to  determine  the 
fluid  velocity  dependence  on  projectile 
speed,  the  bulk  of  the  results  pertain 
to  maximum  projectile  velocity  of 

6.8  ms”1,  achieved  over  a  300  am  travel 
with  the  initial  chamber  pressure  of 

2.8  bars.  The  axial  velocity  traces,  or 
values  which  are  considered  to  be  the 
predominant  component  of  velocity,  were 
obtained  from  different  shots  and  the 
repeatability  of  experiments  was 
generally  within  1%. 

Figure  2  shows  typical  travel-time  and 
velocity-time  curves  of  the  projectile 
obtained  with  as  initial  gas  pressure 
of  2.8  bars.  The  projectile  velocity 
increases  during  the  first  30  ms, 
corresponding  to  about  140  an  travel  of 
the  projectile,  under  the  pressure  forces 
exerted  os  its  base.  The  velocity 
attalns1an  almost  constant  value  of  about 

6.6  ms'1  during  the  next  20  ms  and  there- 
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after  retards  gradually  as  tbe  pressure 
forces  diminish  due  to  expansion  and  tbe 
resistive  forces, due  to  friction,  take 
over.  Tbe  abrupt  fall  In  tbe  projectile 
velocity  observed  after  about  300  on 
travel  58  as)  is  associated  with 
tbe  iaplngeoent  of  tbe  projectile  on  tbe 
compression  spring  wblcb  tben  colls  and 
gradually  retards  the  projectile.  Tbe 
unaffected  length  of  300  mb  travel  will, 
therefore,  be  considered  from  noe  on  as 
tbe  projectile  travel  corresponding, 
for  tbis  particular  initial  pressure,  to 
a  cycle  time  of  56.7  as. 

Tbe  velocity-travel  curves  of  tbe 
projectile  obtained  sltb  various  initial 
pressures  ranging  from  1.7  to  3.B  bars 
are  shown  in  Figure  3.  Increases  in 
initial  pressure  result  in  increased 
projectile  velocities  alaost  la  a 
parabolic  manner  as  oae  would  expect 
from  tbe  energy  balance  which  neglects 
friction.  4s  tbe  sliding  friction 
remains  almost  constant,  tbe  velocity- 
travel  curves  on  tbe  second  half  of  the 
cycle  gradually  ebange  slope  from  a 
negative  value  to  a  positive  (le  from 
deceleration  to  acceleration)  with 
increasing  pressure  forces  as  indicated 
in  Figure  3. 

Tbe  corresponding  traces  of  centreline 
velocity  obtained  inside  the  initial 
volume  and  near  the  projectile  base  at 
t  -  0  are  shown  for  the  gun  cycle  in 
Figure  4.  Tbe  trends  are  similar  to 
those  of  the  projectile  velocity  and 
suggest  that  fluid  velocities  should  scale 
linearly  with  projectile  speed.  Tbis  is 
verified  by  Figure  5  which  Indicates  tbe 
linear  dependence  of  tbe  fluid  velocity 
on  tbe  projectile  speed  obtained  with 
different  initial  gas  pressures.  Is  line 
with  the  Increase  of  overall  fluid 
velocities,  tbe  velocity  fluctuations 
also  Increase  with  increasing  pressure 
and  results  In  similar  turbulence 
intensities  (u/0)  in  each  case  wltb 
maximum  values  of  no  more  than  4%. 

Figure  6  shows  tbe  centreline  velocity 
traces  at  various  axial  locations 
together  with  tbe  projectile  velocity 
obtained  wltb  tbe  initial  pressure  of 
2.6  bars.  Tbe  axial  locations  of  z  •  100 
and  165  mm  correspond  to  the  initial 
chamber  volume  filled  wltb  gas  at  rest. 

The  fluid  velocities  measured  at  these 
locations  therefore  start  from  sero  and 
Increase  in  time  at  a  rate  which 
approachea  tbe  projectile  acceleration 
(do  /dt)  at  locations  clone  to  tbe  base. 

4s  Tbe  projectile  moves  away  from  these 
locations  with  a  continuous  exchange  of 
momentum  between  tbe  fluid  and  tbe 
projectile,  tbe  velocities  decrease  in 


magnitude.  4t  locations  which  are 
initially  occupied  by  tbe  projectile,  tbe 
fluid  has  a  velocity  equal  to  that  of  tbe 
projectile  at  tbe  time  it  crosses  this 
location,  but  tbe  fluid  subsequently  lags 
behind  tbe  faster  moving  projectile  with 
decreasing  velocities  in  tbe  rest  of  tbe 
cycle  as  shown  in  Figure  6.  although  tbe 
maximum  rms  velocity  fluctuations  (5)  in 
tbe  cycle  increases  from  about  Q  08  ms"1 
at  z  •  100  mm  to  about  0.21  ms~  at 
x  -  295  am,  tbe  intensity  remains 
virtually  tbe  same  at  every  location  at  a 
maximum  value  of  about  4%.  It  is  probable 
that  these  fluctuations  stem  mainly  from 
turbulence  production  due  to  normal 
stresses  rather  than  shear  since,  as  will 
be  shown  later,  insignificant  velocity 
gradients  exist  in  most  regions  of  tbe 
flow.  During  tbe  early  part  of  tbe  cycle, 
however,  so  velocity  fluctuations  were 
observed  due  to  tbe  stabilising  effect  of 
strongly  accelerating  flow. 

The  distributions  of  centreline  velocity 
are  shown  in  Figure  7  for  various  times. 
Tbe  most  striking  feature  of  these 
results  is  that  at  any  given  time  tbe 
velocities  are  linearly  distributed  in 
the  space  between  tbe  breech  and  the 
projectile  base  with  the  values  Increasing 
from  zero  at  the  breecb  to  that  of  tbe 
projectile  velocity  at  the  base.  4  direct 
implication  of  this  result  is  that  tbe 
linear  velocity  assumption  widely  used 
by  interior  ball isti clans  for  the 
conversion  of  breecb  pressure  to  base 
pressure,  see  references  1  and  6  for 
example,  is  a  correct  one  at  least  in 
tbe  limiting  case  of  a  single-phase  flow. 

4  description  of  the  flow  field  prevailing 
in  tbe  gun  tube  during  tbe  acceleration 
part  of  tbe  cycle  is  given  in  Figure  B 
which  presents  tbe  radial  profiles  of 
axial  velocity  at  various  times  and  axial 
locations.  Tbe  profiles  reveal  that  the 
Velocities  Increase  wltb  time  and  axial 
location,  and  that  they  are  uniformly 
distributed  across  tbe  cylinder  indicat¬ 
ing  tbat  tbe  flow  is  one-dimensional 
tbrougbout  tbe  space  and  cycle  apart  from 
tbe  thin  regions  near  tbe  wall  boundaries. 

More  detail  of  tbe  flow  sear  tbe  wall 
boundaries  it  i  •  J  (76.7  mm)  and 
2D  (153.4  ms)  is  given  in  Figures  9  and 
10  respectively,  wblcb  presents  tbe 
velocity  traces  at  radial  distances  r 
of  0.5,  1,  2  and  3  mm  away  from  tbe  will. 
These  traces  do  sot  ebon  any  significant 
enhancement  in  velocity  fluctuations  near 
tbe  wall  boundary  and  suggest  tbat  tbe 
boundary  layers  are  laminar  and  remain 
laminar  throughout  tbe  cycle.  Tbe 
velocity  values  at  r  »  0.5  b  are 
generally  of  tbe  order  of  95%  of  tbe 
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velocities  it  r.  •  3  ■  sod  suggest  that 
tile  boundary  lifer  thickness  j  is  less 
than  1  an  corresponding  to  about  1%  of 
the  cylinder  diaoeter.  The  calculations 
of  reference  3  based  os  laminar  flow 
indicate  that  tbe  displacement  thickness 
t. ,  which  is  roughly  30%  of  the  boundary 
lflyer  thickness  (see  for  example 
reference  7),  Is  generally  of  the  order 
of  about  0.3%  of  th<  diaoeter  correspond¬ 
ing  therefore  to  similar  boundary  layer 
thicknesses  to  those  observed  here. 


4  CONCLUSIONS 

Tbe  main  findings  of  this  investigation 

can  be  sunsarlsed  as  follows: 

1.  The  flow  at  any  given  time  in  tbe 
cycle  is  near  one-dimensional  and 
laminar  with  velocity  fluctuations 
corresponding  to  maximum  Intensities 
of  about  4%. 

2.  The  boundary  layers  remain  laminar 
during  the  cycle  and  correspond  to 
about  1%  of  the  cylinder  diameter. 

3.  The  fluid  velocities  are  linearly 
distributed  in  the  space  between 
the  breech  and  projectile  base  with 
velocities  increasing  from  zero  at 
the  breech  to  that  of  the  projectile 
near  the  base. 

4.  The  fluid  velocities  show  a  linear 
dependence  on  the  projectile  speed 
obtained  with  different  initial  gas 
pressures . 
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ABSTRACT 


Time-resolved  and  ensemble  mean  and  rms  fluid  velocities 
obtained  by  laser  Doppler  anemometry  in  a  subsonic  gun  simul¬ 
ator  with  an  inert/single-phase  flow  and  a  projectile  exit 
velocity  of  40m/s  are  reported.  The  results  presented  also 
include  the  time-resolved  measurements  of  projectile  velocity 
and  breech  pressure. 

The  results  show  that  the  calculated  projectile  velocities 
based  on  measured  breech-pressure  record  are  almost  exactly 
equal  to  the  measured  values  and  confirm  that  the  friction  of 
the  system  is  negligible.  Fluid  velocity  traces,  together  with 
ensemble  mean  and  rms  velocity  profiles,  indicate  the  presence 
of  a  wall  boundary  layer  which  occupies  up  to  20%  of  the  tube 
radius  and  with  a  growth  rate  similar  to  steady  turbulent 
boundary  layers.  The  nature  of  velocity  fluctuations  and  the 
shape  of  velocity  profiles  inside  the  boundary  layer  also 
suggest  that  the  boundary  layers  in  the  wake  of  the  projectile 
are  turbulent. 


LIST  OF  SYMBOLS 


A  cross-section  area  of  projectile 

a  local  speed  of  sound 

d  tube  diameter 

friction  force 

g  acceleration  of  gravity 

m  mass  of  projectile 

Pfl  atmospheric  pressure  (back  pressure) 

Pfa  breech  pressure 

PQ  absolute  initial  pressure 

R  tube  radius 

r  distance  from  cylinder  wall 

t  time 

U  ensemble  mean  axial  velocity 

Uc  centreline  mean  velocity 

u'  ensemble  rms  axial  velocity 

Vp  projectile  velocity 

x  Eulerian  coordinate  measured  from  breech  (fluid) 

2  Lagrangian  coordinate  (projectile) 

zc  initial  chamber  length 

Y  isentropic  coefficient 

6  boundary  layer  thickness 

PQ  initial  fluid  density 
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1.  INTRODUCTION 


The  flow  field  behind  the  projectile  in  a  real  gun  is 
characterized  by  high  velocities  of  solid  and  gas  phases  and 
by  exothermic  reactions,  with  temperatures  up  to  3000  K, 
reference  1,  and  unsteady  boundary  layers  at  the  walls  behind 
the  projectile.  The  short  duration  of  the  projectile  travel 
within  the  barrel  adds  to  the  difficulty  of  accurate  measure¬ 
ments  of  the  high  speed  transient  phenomena.  A  related  prog¬ 
ram  of  research  has  recently  been  recommended,  reference  2, 
and  involves  experiments  with  simulators  that  permit  the 
generation  of  well  controlled  flows,  starting  with  more 
simplified  single-phase  inert  flows  and  advancing  progressive¬ 
ly  towards  more  complex  two-phase  reacting  flows.  A  strong 
interest  is  also  being  shown  by  computational  communities 
(references  3 and  4)  which,  it  is  hoped,  will  provide  means  of 
interpolation  and  extrapolation  of  the  experimental  results. 
Heiser  (reference  4),  for  example,  has  obtained  encouraging 
agreement  with  the  experimental  results  of  reference  5,  by 
solving  the  full  Navier-Stokes  equations  for  laminar  flows. 

The  work  presented  here  describes  inert  single-phase 
measurements  of  pressure,  projectile  and  fluid  velocities  in 
a  subsonic  gun  simulator.  The  model  gun  is  of  the  preburned- 
propellant  and  constant-diameter  type, in  ballistic  phrase¬ 
ology.  It  constitutes  a  continuation  and  improvement  of  the 
work  reported  in  reference  5.  The  wall  shear  layers  which 
existed  in  a  gas  expanding  behind  a  projectile  of  reference 
5  were  found  to  be  less  than  1%  of  the  tube  diameter  and  it 
was  difficult  to  provide  detailed  results  which  would  clarify 
the  nature  of  the  unsteady  boundary  layer.  In  the  present 
study  the  length  of  the  travel  of  the  projectile  has  been 
extended  thus  allowing  the  boundary  layers  to  grow  more  than 
previously.  Little  is  known  about  transient  boundary  layers 
(references  6  and  7)  and,  consequently,  particular  attention 
is  paid  to  the  nature  of  the  ensuing  boundary  layers  since 
their  behaviour  accounts  for  momentum  and  heat  transfer  to 
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the  wall  of  the  barrel.  In  addition,  a  novel  optical  tech¬ 
nique  has  been  developed  and  used  successfully  to  obtain  the 
projectile  velocity  as  a  function  of  space  and  time.  This 
work  also  represents  a  stage  in  a  program  of  research  aimed 
at  measuring  in  high-speed  two-phase  flow  configurations  by 
evaluating  measurement  techniques. 


2.  EXPERIMENTAL  TECHNIQUES 
2 . 1  Flow  configuration 

The  flow  arrangement  shown  in  Figure  1  is  similar  to 
that  in  reference  5  but  has  a  projectile  travel  of  I4d  as 
opposed  to  4d  and  a  lighter  projectile  to  achieve  higher 
velocities  with  the  range  of  initial  pressures  previously 
used. 

The  projectile  is  secured  in  position  inside  a  tube  of 
76.5mm  diameter  which  is  pressurized  by  compressed  nitrogen. 
After  the  gas  becomes  quiescent,  the  projectile  is  released 
and  travels  a  distance  of  about  lm  (=13d)  inside  the  tube 
which  is  machined  to  high  precision.  The  initial  chamber  of 
306mm  (4d)  in  length  is  fabricated  from  plexiglass  to  provide 
optical  access  for  the  laser  Doppler  anemometer.  The  remain¬ 
der  of  the  barrel  is  made  of  mild  steel,  with  four  small 
circular  plexiglass  windows,  as  shown  in  Figure  1,  to  allow 
velocity  measurements  at  axial  distances  of  approximately  7, 
10,  13  and  16d  from  the  breech.  A  lock  mechanism  has  been 
designed  to  secure  the  projectile  under  the  pressure  load  at 
its  starting  position  and  to  release  it  smoothly  and  spon¬ 
taneously  without  causing  jitter.  The  initial  sealing  of 
the  chamber  is  achieved  by  means  of  an  0-ring  fitted  around 
the  inner  surface  of  the  tube,  against  which  the  projectile 
rests  at  its  start  position.  After  the  projectile  is 
released  and  during  its  travel  inside  the  barrel  two  plastic 
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rings  fitted  on  the  projectile  prevent  blow-by  of  the  com¬ 
pressed  gas.  The  projectile  is  2  5Cknm  long  and  weighs 
0.953kg.  The  top  and  bottom  parts  are  made  of  aluminium 
alloy  and  a  central  plexiglass  section  houses  a  linear  grati¬ 
cule  for  projectile  velocity  measurements.  At  the  end  of  its 
travel,  the  projectile  is  decelerated  by  compression  of  the 
air  inside  the  lower  part  of  the  barrel.  A  control  valve 
fitted  at  the  bottom  and  connected  to  a  compressed  air  supply 
is  used  to  return  the  projectile  to  its  lock  and  start  pos¬ 
ition  at  the  top  of  the  barrel. 


2 . 2  Measurement  systems 

Emphasis  has  been  placed  on  the  development  of  an 
improved  method  for  measuring  the  projectile  velocity.  The 
chosen  optical  technique  makes  use  of  a  linear  graticule  and 
eliminates  the  need  for  the  alternative  and  mechanically  com¬ 
plex  rack  and  gear  mechanism  of  reference  5  and  removes  the 
ambiguity  in  the  associated  friction  force.  The  graticule 
inside  the  projectile  consists  of  25  lines  per  inch.  The 
measuring  systems  are  shown  in  Figure  1. 

For  the  projectile  velocity  measurements,  the  instrum¬ 
entation  comprises  a  5mW  He-Ne  laser,  a  lens  to  focus  the 
beam  at  its  waist  on  the  graticule,  a  pinhole  and  another 
lens  to  focus  the  emerging  beam  onto  a  photodiode.  The  lens- 
pinhole  combination  reduces  the  field  of  view  of  the  diode 
to  a  minimum  possible  for  the  system  operation  and  thereby 
reduces  the  sensitivity  of  the  system  to  unwanted  optical 
noise  from  reflections  and  flares.  The  graticule  moves  with 
the  projectile  and  transmits  a  pulse  train  whose  period 
coincides  with  the  chopping  of  the  laser  beam  by  the  grati¬ 
cule.  The  output  signal  of  the  diode  is  converted  to  a 
square  wave  by  a  schmitt  trigger,  which  is  competible  with 
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digital  logic.  The  individual  periods  are  timed  by  a  high 
frequency  clock  and  stored  in  a  buffer  memory  from  which 
the  velocity  of  the  projectile  can  be  computed  during  the 
part  of  the  cycle  corresponding  to  the  transient  time  of  the 
projectile  across  the  selected  window.  By  repeating  the 
measurements  at  each  of  the  four  windows,  velocity-travel 
curve  for  the  whole  cycle  can  be  assembled. 

For  the  measurements  of  fluid  velocity  a  dual-beam  laser 
Doppler  anemometer  operated  in  the  forward-scatter  mode  was 
used.  The  anemometer  consisted  of  a  32mW  He-Ne  laser,  a 
focussing  lens,  a  rotating  diffraction  grating,  a  collimating 
lens,  a  collecting  lens  and  a  photomultiplier  as  also  shown 
in  Figure  1.  The  principal  characteristics  of  the  optical 
system  are  given  in  Table  1.  Silicone  oil  particles  were 


TABLE  1  Principal  characteristics  of  optical  system 


Half  angle  of  intersection  of  beams  6.10  degrees 


Fringe  spacing 

Length  of  probe  volume  at  e' 
intensity  level 

Diameter  of  probe  volume  at  e 
intensity  level 

_2 

Number  of  fringes  within  e 
intensity  level 


2.99  urn 


0.908  mm 


0.097  mm 


34 


used  as  scattering  particles  and  were  introduced  in  the  ini¬ 
tial  chamber  prior  to  charging  of  the  volume  at  the  required 
release  pressure.  The  output  of  the  photomultiplier  was 
band-pass  filtered  to  remove  the  pedestal  and  the  high  frequency 
noise  and  input  to  a  purpose  built  frequency  counter  with  an  8/16 
validation  logic  and  a  cycle  time  of  45visec.  The  frequency  counter  was 
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interfaced  to  a  microcomputer  which  sampled  the  digital  out¬ 
put  of  the  counter  and  stored  it  in  memory  at  a  rate  of 
=10KHz.  The  time  information  corresponding  to  each  frequency 
count  was  obtained  from  a  programmable  timer  module  which 
was  activated  to  time  at  a  rate  of  477KHZ  after  the  first 
2mm  of  the  projectile  travel.  The  data  was  subsequently 
processed  by  the  microcomputer  to  determine  and  display  the 
variation  of  velocity  as  a  function  of  time. 

The  initial  pressure  used  was  8  bars  and  achieved  with 
low  grade  nitrogen  taken  from  a  pressurised  cylinder.  The 
initial  pressure  and  the  breech  pressure  record  during  the 
cycle  were  measured  with  a  precision  of  -0.002  bars,  using 
a  piezoelectric  transducer  and  a  charge  amplifier  whose 
analogue  output  was  digitized  and  interfaced  with  the  micro¬ 
computer. 


3.  RESULTS  AND  DISCUSSION 


3 . 1  Projectile  dynamics 

Figure  2  shows  the  measured  velocity-travel  curve  of 
the  projectile  together  with  the  curve  calculated  from  the 
measured  breech  pressure-time  record.  The  measured  traces 
were  found  to  be  repeatable  within  1%  from  cycle  to  cycle 
over  a  large  number  of  cycles  and  the  problem  of  varying 
friction  encountered  with  the  previous  experimental  set-up 
of  reference  5  has  been  removed.  The  projectile  accelerates 
rapidly  from  rest  to  reach  its  exit  velocity  of  about  40m/s 
In  about  30msec.  The  acceleration  of  the  projectile  is 
greatest  at  the  start  of  the  motion  and  decreases  monotonic- 
ally  to  reach  almost  zero  at  the  end  of  the  cycle. 

The  motion  of  the  projectile  can  be  determined  if  the 
pressure  acting  on  its  base  is  known.  The  relation  that 
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gives  the  base  pressure  as  a  function  of  the  breech  pressure 
acting  along  a  characteristic  line  is  (reference  8) : 
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The  projectile  motion  as  a  function  of  the  breech 
pressure  is  then  given  by  the  equations: 
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In  equation  2  the  compressibility  effect  on  the  atmospheric 
side  has  been  ignored  since  under  the  present  conditions  it  does 
not  contribute  more  than  0.5%  to  the  exit  projectile  velocity 
(see  reference  9  for  example).  Preliminary  pressure  measure¬ 
ments  obtained  near  the  muzzle  verified  this  assumption. 

Equations  2  and  3  constitute  a  coupled  set  of  simult¬ 
aneous  ordinary  differential  equations  in  z  and  Vp  in  terms 
of  t  and  can  be  solved  numerically,  by  a  standard  fourth 
order  Runge-Kutta  Algorithm.  The  friction  force  term  in 
Equation  2  is  taken  to  be  equal  to  zero  since  no  a  priori 
model  for  it  exists.  The  calculated  values  for  Vp  in  Figure 
2  indicate  a  value  at  the  exit  position  of  40.0m/s  which  is 
almost  exactly  equal  to  the  measured  value  and  thus  suggests 
a  negligible  friction  force.  This  fact  assumes  particular 
significance  for  numerical  predictions,  which  obviates  the 
need  for  considerations  of  the  wall  friction  force. 


jk 
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3. 2  Fluid  velocity  measurements 

The  centreline  velocity  traces  of  the  fluid  together 
with  the  projectile  velocity  measurements  are  shown  in  Figure 
3  for  the  axial  stations  of  x  =  3d,  x  =  7.27d,  x  =  10.20d, 
x  =  13.13d  and  x  =  16.07d.  At  x  =  3d,  which  corresponds  to 
the  position  inside  the  initial  chamber,  the  fluid  velocity 
increases  monotonically  with  time  to  reach  a  maximum  of 
14m/s  at  around  10msec.  During  the  acceleration  period  no 
fluctuations  are  present.  The  velocity,  after  reaching  the 
maximum,  decreases  to  7m/s  on  average  at  about  35msec  and 
fluctuations,  though  of  low  magnitude,  start  to  appear  with 
decreasing  acceleration.  The  fluid  velocities  immediately 
behind  the  projectile  have  been  measured  at  subsequent  axial 
stations  and  increase  in  accord  with  the  projectile  velocity 
which  forms  an  envelope  to  the  fluid  velocities.  At  these 
axial  stations,  the  fluid  velocities  decay  in  an  almost 
linear  fashion  with  time,  as  the  projectile  moves  away.  The 
fluctuations  increase  with  axial  distance. 

Figure  4  shows  fluid  velocity  traces  at  x  =  7.27d  for 
different  radial  positions  from  the  wall  and  suggests  that 
fluid  is  turbulent  and  certainly  far  from  laminar.  All 
radial  positions  reveal  similar  decay  in  almost  linear 
fashion  with  time.  The  fluid  velocities  at  13msec  are  almost 
the  same  as  the  corresponding  projectile  velocity  of  30m/s  as 
expected.  As  the  wall  is  approached,  the  relative  magnitude 
of  the  fluctuations  increases  and  fewer  velocity  values  are 
recorded.  This  suggests  that  the  density  of  scattering 
particles  is  lower  than  in  the  core  of  the  flow.  It  can 
also  be  seen  that  the  mean  level  of  the  instantaneous  vel¬ 
ocities  drops  close  to  the  wall. 

In  order  to  build  a  picture  of  the  flow  inside  the 
barrel  at  a  certain  time,  ensemble  mean  and  rms  values  need 
to  be  considered  during  a  finite  time  interval  and  over  many 
single  realizations.  The  experimental  sample  size  was  chosen 
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to  be  greater  than  80  within  a  time  interval  of  3msec  as  a 
compromise  between  the  need  to  reduce  the  time  gradient 
broadening  and  the  number  of  samples  necessary  to  obtain  an 
acceptable  statistical  error.  The  number  of  at  least  80 
samples  ensured  a  maximum  statistical  error,  according  to 
reference  10, in  the  mean  velocity  of  2%.  In  the  core  of 
the  flow  where  the  sample  sizes  were  generally  greater  than 
120,  the  error  is  less  than  1%.  The  corresponding  error  in 
the  rms  velocity  is  of  the  order  of  10%. 

Other  sources  of  error  are  gradient  broadenings  due  to 
the  variation  of  the  velocity  within  the  averaging  time 
interval  of  3msec  and  within  the  measuring  probe  volume. 

The  maximum  possible  error  including  the  accuracy  of  the 
measurement  system,  gradient  broadenings  and  the  statistical 
uncertainty  is  about  4%  in  the  mean  velocity  and  less  than 
20%  in  the  rms  velocity. 

Figure  5  presents  radial  profiles  of  the  mean  and  rms 
axial  velocity  inside  the  barrel  for  a  time  of  26msec  at 
x  *  7.27d,  x  =  10.20d  and  x  =  13.13d  and  represents  the 
state  of  the  fluid  motion  when  the  projectile  has  almost 
reached  the  muzzle.  The  mean  velocity  profiles  show  similar 
trends,  with  a  flat  core  which  occupies  almost  80%  of  the 
barrel  radius.  Near  the  wall  both  mean  and  rms  profiles 
reveal  the  existence  of  a  shear  layer.  The  rms  values  are 
low  over  the  bulk  of  the  flow,  around  4%  of  the  local  mean 
velocity,  and  increase  to  reach  10%  of  the  mean  velocity 
near  the  wall.  Also  shown  in  Figure  5  is  the  axial  distrib¬ 
ution  of  the  boundary  layer  thickness,  which  increases  from 
0. 18R  at  x  *  7.2"'d  to  0.21R  at  x  *  10.20d  and  decreases  to 
O.ISR  at  x  ■  13.13d  and  should  go  to  zero  at  the  projectile 
base. 

Little  is  known  about  the  behaviour  of  transient  boun¬ 
dary  layers  and  a  comparison  with  steady  boundary  layers  may 
be  helpful.  Figure  6a  shows  the  variation  of  the  measured 
boundary  layer  thickness  with  axial  distance  at  a  fixed  time 
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of  26  msec  together  with  the  calculated  variations  of  a  steady 
incompressible  laminar  boundary  layer  starting  at  the  breech 
and  steady  incompressible  turbulent  boundary  layers  with  zero 
and  favourable  pressure  gradient  (reference  11) ,  starting  at 
the  end  of  the  initial  chamber  where  the  boundary  layers  are 
believed  to  be  laminar,  reference  5.  For  the  calculations  of 
boundary  layer  thickness  with  favourable  pressure  gradient,  a 
linear  variation  of  core  velocity  between  the  breech  and 
projectile  base  was  used  at  the  chosen  time  to  evaluate  the 
corresponding  pressure  gradient.  The  measured  boundary  layer 
thicknesses  are  more  closely  related  to  the  turbulent  boundary 
layers  as  shown  in  the  figure. 

Figure  6b  shows  the  measured  velocity  profiles  inside 
the  boundary  layer  for  different  times  and  axial  locations, 
together  with  1/7-power  law  of  turbulent  steady  flow  inside 
a  tube.  The  measured  velocity  prefiles  follow  closely  the 
1/7-power  law.  The  exponent  for  the  measured  profiles  can 
be  found  by  the  least-squares  technique  which  give  values 
from  1/6.6  for  the  lowest  bulk  Reynolds  number  based  on  the 
diameter  of  the  tube,  which  occurs  at  the  earliest  parts  of 
the  projectile  motion  to  1/9.9  for  the  highest  bulk  Reynolds 
number  which  occurs  when  the  projectile  reaches  the  muzzle. 

The  calculated  shape  factor  ranged  from  1.2  to  1.3  and  again 
suggests  that  the  measured  boundary  layers  are  turbulent 
(reference  12) . 

Figure  7  shows  the  profiles  of  the  mean  axial  velocity 
at  x  =  7.27d  for  three  different  times  during  the  cycle  of 
14.5,  26  and  30.5msec  and  represents  the  evolution  of  the 
mean  axial  flow  with  time  at  a  fixed  position  (Eulerian 
frame) .  It  can  be  seen  that  the  boundary  layer  thickness 
increases  from  0.12R  at  14.5msec  to  reach  0.19R  later  in  the 
cycle.  The  boundary  layer  at  a  fixed  position,  therefore, 
grows  with  time  during  the  cycle  and  reaches  an  almost  con¬ 
stant  thickness.  Figure  8  shows  the  evolution  of  the  mean 
and  rms  axial  velocity  profiles  at  a  fixed  location  of  60mm 
behind  the  moving  projectile,  so  that  it  represents  the 
flow  following  the  motion  of  the  projectile  (Lagrangian 
frame).  In  common  with  Figure  7,  the  profiles  at  14.5msec 
show  a  boundary  layer  thickness  of  0.12R  which  grows,  however, 
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at  a  much  slower  rate  to  reach  0.13R  at  20.5msec  and  0.15R 
at  26msec. 


4 .  CONCLUSIONS 


The  following  main  conclusions  may  be  extracted  from 
the  preceding  text: 

1.  The  relationship  between  the  breech  pressure  and 
the  projectile  velocity  may  be  determined  by  consideration 
of  one-dimensional,  unsteady  forms  of  the  continuity  and 
momentum  equations.  The  measured  breech  pressure  record 
leads  to  calculated  projectile  velocity  values  within  0.5% 
of  the  measured  ones  and  confirms  that  the  friction  between 
the  projectile  and  the  tube  wall  is  negligible. 

2.  The  fluid  velocity  traces,  together  with  the  corres¬ 
ponding  ensemble  mean  and  rms  velocity  profiles,  indicate  the 
presence  of  a  shear  layer  which  exists  behind  the  projectile 
and  occupies  around  20%  of  the  tube  radius.  In  the  initial 
chamber,  the  boundary  layer  thickness  is  of  the  order  of  1% 
of  the  radius. 

3.  The  nature  of  the  fluid  velocity  fluctuations  and 
the  shape  of  the  ensemble  mean  and  rms  velocity  profiles  in¬ 
side  the  boundary  layer,  suggest  that  the  boundary  layer  is 
turbulent.  The  evolution  of  the  boundary  layer  thickness  as 
compared  to  steady  turbulent  and  laminar  boundary  layers 
shows  that  it  is  closer  to  the  former. 
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Figure  Is  Schematic  diagram  of  experiJtental  set  up 


Figure  4:  Fluid  velocity  traces  at  x  ■  7.27d  for  different  radial 
positions  frcm  the  veil 


d/9 


0 


10 


30 


20 

U  &  u'  in  ml s 


in 


t=26  msec 


Figure  5:  Mean  and  rms  velocity  profiles  at  t  *  26msec 
x  =  7.27d,  x  «  10.20d  and  x  *  13.13d  . 


t  =  26  ms 


r/R  r/R  r/R 


Figure  7:  Mean  velocity  profiles  at  x  =  7.27d  for  different  times 


Figure  8:  Mean  and  rms  velocity  profiles  at  x  =  z  +  zo  - 
for  different  times 
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Abstract 


Fluid  and  solid  particles  (500  im  diameter)  velocities  obtained  by 
laser  Doppler  anemooetry  in  a  subsonic  gun  simulator  with  inert  single-  and 
two-phase  flow  are  reported.  The  gun  simulator  is  of  the  preburned 
propellant  chambered  type  with  a  gradual  contraction  in  the  top  initial 
chamber  which  incorporated  a  particle  injector  mechanism.  The  results  show 
that  the  effect  of  the  chambered  geometry  on  fluid  and  projectile  velocity 
and  on  the  growth  of  the  wall  shear  layer  remained  negligible  within  the 
subsonic  range  considered.  Comparison  between  the  velocity  of  500  iaa 
DiaXon  particles  and  that  of  the  fluid  obtained  separately  in  the  absence 
of  solid  particles  revealed  that  the  particles  lag  behind  the  fluid  in 
regions  closer  to  the  projecile  because  of  their  relatively  large  inertia, 
and  that  the  lag  increases  with  axial  distance  from  the  breech  reaching  a 
measured  maximum  of  about  2%  near  the  end  of  the  cycle.  The  observed  rms 
levels  of  particle  velocity  are  higher  than  those  of  the  fluid  in  the  bulk 
of  the  flow  but  attain  similar  values  near  the  tube  wall. 
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In  a  real  gun  environment,  the  interior  ballistic  flow  ia 
transient  and  reactive  with  fast  burning  propellant  grains,  metal  and 
eoot  particles  of  various  sire*  and  of  irregular  shapes.  The  two-phase 
nature  of  the  flow  adds  to  its  complexity  and  makes  it  almost 
impossible  to  probe  the  flow  for  velocity  information.  Progress  in 
modelling  and  calculation  methods  (1,  2),  on  the  other  hand,  calls  for 
understanding  of  the  Characteristics  of  unsteady,  two-phase  flows  in 
relevant  arrangements,  and  the  corresponding  measurement  problems  imply 
that  simplified  geometries  are  required. 

The  flowfteid  behind  a  moving  projectile  in  a  single-phase 
subsonic  gun  simulator  has  been  quantified  (3,  4)  and  the  related 
transient  boundary  layers  were  shown  to  be  similar  to  steady  turbulent 
boundary  layers.  This  paper  reports  experiments  with  inert  single-  and 
two-phase  flows  in  a  similar  gun  simulator  and  represents  a  further 
step  towards  reality.  The  single-phase  flow  experiments  quantify  the 
effect  of  initial  chamber  geometry  including  a  gradual  contraction 
similar  to  those  of  real  gun  systems  and  the  two-phase  flow  results 
quantify  the  behavior  of  500  fm  Diakon  particles.  The  latter  make  use 
of  a  particle-injector  mechanism  fitted  to  the  breech  Chamber  to 
introduce  the  annodispersed  particles  with  initial  volume  concentration 
of  about  1.3%. 

The  following  section  describes  the  gun  simulator  and  the 
experimental  techniques.  Section  3  presents  and  discusses  the  results 
and  Section  4  summarises  the  main  findings  of  this  wsrk. 

Experimental  System 

A  diagram  of  the  gun  simulator,  which  is  mounted  vertically,  is 
shown  in  Figure  1  and  is  similar  to  that  of  (3),  with  the  exception  of 
the  initial  Chamber  geometry  tdiiCh  ooeprlses  a  gradual  contraction  from 
a  diameter  of  150  to  76.5  am.  The  initial  chamber  of  (3)  had  a 
constant  diameter  of  76.5  am  and  approximately  the  same  volume .  The 
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present  arrangement  simulates  sore  closely  the  Initial  chamber  geometry 
of  real  gun  systems  (5)  and  incorporated  a  13  m  thick  disk,  with 
twelve  14  ■e-diaaeter  holes  distributed  as  shewn  on  figure  1.  and 
initially  filled  with  500  tm  Dlakon  particles.  The  holes  were  sealed 
by  spring-loaded  cones  pushed  against  them  by  the  projectile  in  its 
lock  position.  On  release  of  the  projectile,  the  oones  retracted  and 
the  particles  entered  the  initial  chasber.  The  number  of  holes,  their 
distribution,  cone  angle  and  the  particular  initial  chamber  geometry 
were  Chosen  to  cause  the  particle  distribution  in  the  Initial  volume  to 
be  as  uniform  as  possible.  n>e  volume  above  the  disk  was  less  than  5% 
of  the  total  initial  volume  with  consequently  small  jet  momentum.  With 
this  arrangement  an  initial  particle  concentration  of  up  to  1.3%  by 
volume  was  achieved  and  corresponded  to  a  mass  loading  of  about  2.3. 

The  projectile  was  secured  in  position  at  the  top  of  the  barrel  by 
a  lock/release  mechanism  as  shewn  in  Figure  1,  and  the  initial  chamber 
pressurised  by  compressed  nitrogen.  after  the  gas  became  quiescent, 
the  projectile  was  released  and  the  particles  injected.  The  projectile 
travelled  a  distance  of  about  1  m  inside  the  76.5  am  diameter  tube 
which  was  machined  to  high  precision  and  contained  four  small 
plexiglass  windows  to  provide  optical  access  for  laser  Doppler 
anousmetry  measurements  at  axial  locations  of  approximately  5d,  id,  lid 
and  14d  from  the  breech.  The  initial  sealing  of  the  chamber  was 
achieved  by  means  of  an  O-ring  fitted  around  the  inner  surface  of  the 
tube  against  which  the  projectile  rested  at  its  start  position,  hfter 
the  projectile  was  released  and  during  its  travel  inside  the  tube,  two 
PTFE-rings  fitted  to  the  projectile  prevented  blcw-by  of  the  compressed 
gas.  The  projectile  was  250  ms  in  length  and  weighed  0.953  kg.  The 
top  and  bottom  parts  were  made  from  aluminium  alloy  and  the  central 
plexiglass  section  housed  a  linear  graticule  with  25  lines  per  inch  for 
projectile  velocity  measurements.  At  the  end  of  its  travel,  the 
projectile  was  decelerated  by  the  compression  of  air  inside  the  lower 
part  of  the  barrel.  A  oontrol  valve  fitted  at  the  bottom  and  connected 
to  a  compressed  air  supply  was  used  to  return  the  projectile  to  its 
lock  and  start  position  at  the  top  of  the  barrel. 
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The  projectile  velocity  wee  Measured  with  the  optical  ayetee  of 
Figure  2  and  ooeprieed  a  5  m  He  We  laser,  a  lens  to  focus  the  beaa  at 
its  waist  on  the  graticule,  a  pinhole  and  a  second  lens  to  focus  the 
asmrging  beaa  onto  a  photodiode.  The  graticule  Moved  with  the 
projectile  and  transaitted  a  pulse  train  tdtose  period  coincided  with 
the  chopping  of  the  laser  beaa  by  the  graticule.  the  individual 
periods  were  tiaed  by  a  high  frequency  clock  and  stored  in  a  buffer 
aonory  and  the  projectile  velocity  coaputed  during  the  part  of  the 
cycle  corresponding  to  the  transit  tine  of  the  projectile  across  the 
selected  window.  The  weasureaenta  were  repeated  at  each  of  the  four 
windows  to  provide  a  velocity-travel  curve  for  the  Whole  cycle. 

The  aeasureaents  of  fluid  and  particle  velocity  were  obtained  with 
the  dual-bean,  forward-scatter  laser  Doppler  anemometer  shown  in  Figure 
2.  Silicone  oil  droplets  (mips)  were  used  to  represent  the  gas  flow 
and  the  anoaooeter  was  arranged,  aa  in  (3),  with  a  32  wU  Be-He  laser,  a 
focussing  lens,  s  rotating  diffraction  grating,  a  collimating  lens,  a 
collecting  lens  and  a  photomultiplier .  For  the  measurement*  of 
particle  velocity,  the  anemometer  was  similar  to  that  described  above 
with  a  narrow  slit  in  front  of  the  collecting  lens.  The  Doppler 
Signals  were  processed  by  a  purpose-built  frequency  counter  with  a 
sampling  frequency  of  22  KHz.  The  frequency  counter  was  interfaced  to 
a  microcomputer  which  sampled  the  digital  output  of  the  counter  and 
stored  it  in  memory  at  a  rate  of  about  10  KHz.  The  time  information 
oor responding  to  each  frequency  count  was  obtained  from  a  programmable 
timer  module.  The  data  was  subsequently  processed  by  the  microcomputer 
to  determine  and  display  the  variation  of  velocity  as  a  function  of 
time. 

The  Initial  pressure  was  •  bars  ( gauge )  and  achieved  with  nitrogen 
taken  from  a  preesurisad  cylinder.  The  Initial  pressure  end  the 
breach -pressure  time  record  during  the  cycle  were  measured  with  a 
piezo-electric  transducer  and  a  Charge  amplifier  whose  output  was 
digitised  and  interfaced  to  the  microprocessor. 

A  detailed  analysis  of  the  errors  associated  with  the  experimental 
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technique*  ho*  been  provided  in  (3).  The  aaxisaas  error*  ere  0.5%,  i%, 
4%  and  15%  for  pressure,  projectile  velocity  and  at  an  and  ras  values  of 
fluid  and  particle  velocity,  respectively  and  are  considered  not  to 
affect  the  conclusions. 

3.  Results  and  Discussion 


3.1  Projectile  dvnasu.es 

The  measured  projectile  velocity-travel  curves  obtained  with 
single-phase  flow  and  with  various  initial  chamber  geometries  are  shown 
in  Figure  3.  The  geometries  include  the  gradual  contraction  of  Figure 
1,  the  gradual  contraction  with  the  particle-injector  mechanical  end  the 
oonstant-dlameter  arrangement  of  (3).  the  velocity-travel  curve 
obtained  in  the  presence  of  particles  is  not  included  in  this  figure 
since  it  is  Identical  to  the  single-phase  curve.  Figure  4  shows  the 
corresponding  breech-pressure  time  records  and  that  the  pressure 
variations  of  the  constant  disaster  geometry  are  almost  equal  to  the 
chambered  geometry  with  the  Injector.  The  muscle  velocities  are  45  V* 
with  the  gradual  contraction,  42  a/s  with  the  particle  injector 
Installed  and  40  w/u  with  the  constant-diameter  Chamber  j  these 
differences  stem  mainly  from  the  Changes  in  loading  parameter  (mass  of 
gas/aass  of  projectile )  (6)  associated  with  the  differences  in  the 
Initial  volume  of  the  three  configurations.  The  initial  volume  with 
the  gradual  contraction,  is  14%  larger  than  that  when  the  injector  is 
fitted  and  there  is  22%  difference  between  the  gradual  contraction  and 
the  constant -diameter  volumes.  Calculations  performed  with  the 
one-dimensional  methods  described  in  (3)  using  the  isentropic  variation 
of  static  pressure  at  the  breech  tout  taking  into  account  the 
extra-volume  due  to  chaaberage,  yielded  results  within  2%  of  the 
aeasureaents i  this  confirms  that  the  Initial  geometry  has  negligible 
effect  on  the  projectile  dynamics  for  the  subsonic  flows  considered 
here. 


3.2  Single-phase  flow 


7 


Xn  order  to  quantify  the  effect  of  gradual  contraction  on  flow 
Characteristics,  and  in  particular  on  the  growth  of  boundary  layers, 
the  aean  and  ms  values  of  the  axial  velocity  sere  obtained  at  the 
first  window,  computed  for  a  time  of  14.5  as  and  plotted  in  profile 
form  in  Figure  5,  where  r  denotes  the  radial  position  froei  the  tube 
wall  and  R  is  the  tube  radius.  Wie  figure  also  shows  the  corresponding 
profiles  obtained  with  the  constant-diameter  Chamber  of  ( 3 ) .  With  the 
gradual  contraction,  the  mean  velocities  are  higher  by  about  2  a/a  (7%) 
due  mainly  to  the  difference  in  loading  parameter  and  this  is  also 
evident  in  Figure  6  which  shows  the  axial  distributions  of  centreline 
mean  velocity  for  the  two  chamber  geometries  and  26  me  after  release  of 
the  projectile,  the  fluid  velocities  with  the  gradual  contraction  are 
higher  at  the  axial  locations  corresponding  to  the  first,  second  and 
third  window,  and  the  virtual  origin  established  by  extrapolating  the 
measured  values  to  a  zero  velocity  is  about  -0.92  x/d  (-70  am)  Which, 
corresponds  to  the  extra  volume  associated  with  the  gradual 
contraction. 

The  effect  of  initial  chamber  geometry  on  the  mean  flow 
characteristics  In  the  wake  of  the  projectile  appears  to  be 
Insignificant  and  the  distributions  of  velocity  In  the  initial  volume 
and  in  the  near  upstream  and  downstream  regions  of  the  contraction  are 
far  from  linear.  The  rma  levels  of  the  velocity  fluctuations,  are 
almost  identical  for  the  two  geometries,  with  values  of  around  3%  in 
the  bulk  of  the  flow  and  up  to  9%  in  the  boundary  layer.  The 
corresponding  boundary  layer  thicknesses  for  ths  two  geometries  also 
oompare  well  and  were  found  to  be  12%  and  11%  of  the  radius  of  the  bore 
for  14.5  and  26  ms,  reapectively. 

Fluid  velocities  were  also  measured  with  the  gradual  contraction 
and  in  the  presence  of  the  particle- injector  mechanism  and  Figure  7 
shows  the  profiles  of  mean  and  rma  velocities  obtained  at  the  second 
window  and  at  a  time  of  22  asi  the  profiles  obtained  without  the 
injector  are  also  shown  to  allow  comparisons.  The  profiles  are  almost 
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identical  though  within  the  boundary  layer  the  Injector  caee  displays 
lMer  velocities,  <d\ich  la  surprising  since  the  loading  parameters,  and 
thus  the  projectile  velocities,  are  different.  The  result  can  be 
explained  for  the  core  flow  with  reference  to  Figure  ■  Which  shews  the 
axial  distributions  of  oentreline  eaan  velocity  for  the  two  cases  and 
for  a  ties  of  26  ms .  The  origin  of  the  axial  distance  in  this  figure, 
however,  has  been  chosen  to  correspond  to  the  virtual  origin,  (6),  of 
the  two  configurations,  rather  than  the  geometric  breech-position  of 
Figure  6.  It  is  clear  that  although  the  axial  positions  of  the  windows 
free  the  virtual  origin  (x‘ )  are  different,  the  fluid  velocities  are 
within  5%  of  each  other.  This  result  also  suggests  that  the  injector 
mechanism  does  not  affect  the  as  an  and  turbulent  nature  of  the  flow 
fields  in  the  barrel  and  in  the  wilt  of  the  projectile. 

The  nature  of  the  boundary  layers  with  the  constant  diameter 
Chamber  was  addressed  in  (3)  and  they  were  found  to  be  turbulent  with 
characteristics  similar  to  steady  turbulent  boundary  layers.  The 
boundary  layers  measured  with  the  gradual  contraction  have,  as 
expected,  very  similar  features. 


3.3  Two-c 


se  flow 


The  velocities  of  the  500  jim  particles  were  measured  and  can  be 
compared  with  fluid  velocities  measured  in  the  absence  of  the 
particles.  The  volume  concentration  of  particles  at  the  beginning  of 
the  cycle  was  around  1.3%  and  reduced  to  0.3%  at  the  end  of  the  cycle, 
so  that  the  particle-to-particle  interaction  can  be  assumed  to  be  ssmll 
throughout  the  cycle.  The  results  are  presented  as  ensemble  mean  and 
rms  velocity  profiles  in  Figure  9,  together  with  the  corresponding 
single-phase  flow  profiles.  They  were  obtained  in  the  first,  second 
and  third  windows,  corresponding  to  axial  locations  from  the  breeCh  of 
4.96d,  7.90d  and  10.«4d  respectively,  and  are  presented  at  different 
times  In  the  cycle  due  to  the  variation  in  the  occurrence  and  duration 
of  the  particle  signal  at  different  locations*  in  tbs  tuba.  Sines  the 
fluid  velocities  were  measured  In  the  absence  of  particles,  conclusions 
about  slip  velocity  need  to  be  drawn  with  care.  It  can  be  seen  that 
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the  solid  particles  lag  behind  the  fluid  by  times  which  increase  with 
axial  distance  as  the  projectile  is  approached.  The  difference  between 
the  fluid  and  particle  velocity  increases  in  the  bulk  of  the  flow  free 
about  1  m/»  ( 3% )  at  x/d  -  4.96  and  t  -  20.5  ee  to  9  m/u  (21%)  at  - 
10.94  and  t  “  27.5  as.  This  lag  is  due  to  the  Inertial  effects 
associated  with  solid  particles  in  a  strongly  accelerating  flovf  laid 
particularly  in  regions  close  to  the  projectile  base.  In  this 
Lagrangian  description  therefore,  the  particles  are  always  seen  to  lag 
behind  the  fluid  phase,  in  the  isMdiate  vicinity  of  the  projectile. 
The  particle  res  velocities  are  generally  higher  in  the  bulk  of  the 
flow  than  the  corresponding  single-phase  values  and  the  turbulence 
intensity  at  x/d  »  4.96,  for  exaaple,  is  almost  7%  in  the  bulk  flow  in 
contrast  to  the  single-phase  value  of  about  4%.  the  higher  res  levels 
observed  with  particles  should  be  treated  with  caution,  since  the  fluid 
m  values  were  Measured  in  the  absence  of  particles  and  the  carrier 
phase  should  have  a  different  turbulent  field.  Depending  on  the  else, 
concentration  and  relaxation  tlms  of  the  particles,  the  turbulence 
structure  can  be  Modified  by  the  particles  and  this  nay  Manifest  itself 
as  an  increase  in  turbulence  for  particles  with  relatively  large  else 
and  relaxation  tine  (8  and  9).  It  is  also  probable  that  the  increased 
levels  can  be  associated  with  the  initial  boundary  conditions  ieposed 
in  the  present  experiments .  Although  the  particle-injector  Mechanise 
was  designed  to  alniaise  disturbance  in  the  flow  prevailing  down  the 
tube,  and  this  was  proven  to  be  the  case  for  the  single-phase  flow  as 
previously  discussed,  it  is  possible  that  a  disturbance  superimposed  on 
particle  Motion  will  survive  throughout  the  cycle,  contrary  to  its 
supression  for  the  carrier  phase  due  to  the  stronger  acceleration 
resulting  from  the  expansion  procsss  and  the  gradual  contraction 
Introduced  downs  t  re  as  of  the  injector.  The  particle  res  levels  are 
higher  than  the  fluid  phase  in  the  bulk  of  the  flow,  they  are  very 
sis  liar  to  those  of  the  fluid  near  tha  wall  and  the  effect  of  particles 
on  wall  boundary  layer  thickness  seems  to  be  weak. 

The  Mean  velocity  profiles  of  figure  9  suggest  that  the  particle 
boundary  layers  are  of  the  sans  order  of  thickness  as  ths  fluid  phase 
boundary  layers  have  been  measured  in  the  absence  of  particles. 
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The  variation  in  tha  occurrence  and  duration  of  the  particle  signals  at 
different  location*  In  the  barrel  during  the  cycle,  has  Bade  it 
Impossible  to  obtain  a  oaaplet*  profile  within  the  boundary  layer,  but 
the  measured  turbulence  intensities  are  consistent  with  those  of  the 
single-phase  flow. 

Figure  10  shows  the  oent reline  aean  variation  of  particle  and 
fluid  velocity  with  time  at  x/i  *■  4.96,  in  a  Eulerian  frame.  These 
aean  values  were  evaluated  fix*  around  12  velocity  traces  two  of  «dtiCh 
are  also  shewn  in  Figure  10  for  each  phase.  The  fluid  velocity  is 
higher  at  the  beginning  of  the  cycle  by  about  2  a/e  and  by  about  0.8 
8/*  at  t  ■  23  as.  Later  in  the  cycle  as  the  projectile  moves  away  the 
slip  velocity  Changes  sign,  as  the  particles  catch  tip  and  even  lead  the 
fluid  phase.  The  flowfield  at  this  fixed  location  decelerates  with 
time ,  as  seen  on  Figure  10,  and  the  aean  Reynolds  number  is  also 
continuously  decreasing  with  time.  This  situation  la  in  agreement  with 
the  steady  downflow  of  a  solid/fluid  suspension  where  it  is  knewn  that, 
as  the  Reynolds  number  decreases,  the  lead  of  the  particles  over  the 
fluid  becomes  progressively  larger  (7). 

4.  Conclusions 

The  main  conclusions  of  this  investigation  follow 

(1)  The  effect  of  initial  chamber  geometry  on  projectile  motion  as  well  as 
on  flaw  characteristics  in  the  wake  of  the  projectile  is  found  to  be 
insignificant  within  the  subsonic  range  considered  and  the  relevant 
parameter  is  the  ratio  of  the  mass  of  gas/Bass  of  projectile. 


(2)  The  comparison  between  the  velocity  of  500  i»  Diakon  particles  and  that 
of  the  fluid  obtained  separately  in  the  absence  of  solid  particles 
revealed  that  the  particles  lag  behind  the  fluid  in  regions  closer  to 
the  projectile  because  of  their  relatively  large  inertia  and  that  the 
lag  increases  with  axial  distance  from  the  breech  reaching  a  measured 
maximum  of  ab>ut  20%  near  the  end  of  the  cycle. 
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(3)  ft*  obNtnd  ns  levels  of  particle  velocity  ere  higher  than  thoaa  of 
ths  fluid  in  ths  bulk  of  ths  flow  but  attain  similar  Talus*  naar  the 
tubs  wall.  this  obssrvation  Should,  however,  be  treated  with  care 
since  ths  fluid  as  loci  ties  were  measured  in  ths  absence  of  ths  solid 
particles . 

(4)  at  fixed  locations  in  the  barrel,  the  flowfield  experiences 
deceleration  and,  for  the  saw  reason  as  in  (2)  above,  the  particle 
initially  lags  behind  but  later  in  the  cycle  leads  the  fluid  slightly. 
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figure  2  Schematic  of  experimental  aye tern. 
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Projectile  velocity-travel  curves. 
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Breech-pressure  tine  records. 
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Effect  of  initial  chaster  geometry  on  single-phase  flow 
characteristics . 
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Figure  10 


Centreline  Been  velocities  of  particles  end  fluid  during 
ballistic  cycle  end  it  i  ■  4<96d. 
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